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7^; ABSTRACT 

O . The Euclid mission is the second M-class mission of the ESA Cosmic Vision programme, with 

|H ' the principal science goal of studying dark energy through observations of weak lensing and 

t/^ I baryon acoustic oscillations. Euclid is also expected to undertake additional Legacy Science 

^ programmes. One such proposal is the Exoplanet Euclid Legacy Survey (ExELS) which will 

be the first survey able to measure the abundance of exoplanets down to Earth mass for host 
ff^ \ separations from ^1 AU out to the free-floating (unbound) regime. The cold and free-floating 

^> i exoplanet regimes represent a crucial discovery space for testing planet formation theories. 

\D ■ ExELS will use the gravitational microlensing technique and will detect 1000 microlensing 

events per month over 1 .6 deg^ of the Galactic bulge. We assess how many of these events will 
have detectable planetary signatures using a detailed multi-wavelength microlensing simula- 
tor — the Manchester-Besan§on microLensing Simulator (MAB/^LS) — which incorporates 
\^ , the Besangon Galactic model with 3D extinction. MAB/iLS is the first theoretical simulation of 

' microlensing to treat the effects of point spread function (PSF) blending self-consistently with 

the underlying Galactic model. We use MAB/iLS, together with current numerical models for 
the Euclid PSFs, to explore a number of designs and de-scope options for ExELS, including 
the exoplanet yield as a function of filter choice and slewing time, and the effect of system- 
atic photometry errors. Using conservative extrapolations of current empirical exoplanet mass 
functions determined from ground-based microlensing and radial velocity surveys, ExELS 
can expect to detect a few hundred cold exoplanets around mainly G, K and M-type stellar 
hosts, including ^45 Earth-mass planets and ^6 Mars-mass planets for an observing pro- 
gramme totalling 10 months. ExELS will be capable of measuring the cold exoplanet mass 
function down to Earth mass or below, with orbital separations ranging from ~1 AU out to 
infinity (i.e. the free-floating regime). Recent ground-based microlensing measurements indi- 
cate a significant population of free-floating Jupiters, in which case ExELS will detect hun- 
dreds of free-floating planets. ExELS will also be sensitive to hot exoplanets and sub-stellar 
companions through their transit signatures and this is explored in a companion paper. 

Key words: gravitational lensing: micro — planetary systems — planets and satellites: de- 
tection — Galaxy: bulge — stars: low-mass 

1 INTRODUCTION 
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The discovery of exoplanetary systems is accelerating rapidly with 
over 860 exoplanets confirmed from ground- and space-based ob- 
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servationj3™d another ^2700 cand idates detected with the Kepler 
space telescope jBatalha et alj|2013h . This is providing a wealth of 
knowledge on the distribution function of, primarily, hot exoplanets 
at host separations < 1 AU around FGK-type stars. Recent obser- 
vations by the Kepler space-based transit mission indicate that low- 
mass exoplanets appear to be common and tha t around 20% of stel - 
lar hosts have multiple planets orbiting them (iBatalha et al.ll2013h . 
Resul ts from 8 years of o bservations by the HARPS radial velocity 
team jMavor et alj|201 ll) indicate that half of Solar-type stars host 
planets with orbital periods below 100 days. The frequency of ex- 
oplanets in the Super-Earth to Neptune (SEN) mass range shows 
a sharp increase with declining mass and no preference for host 
star metallicity. HARPS also finds that most SEN planets belong to 
multiple exo planet systems. An analysis by HARPS of its M dwarf 
star sample jBonflls et alj|2013[) indicates that low-mass exoplan- 
ets are also common around low-mass stars and that the fraction rj 

of M dwarf host stars with habitable planets is remarkably high at 

„ _ 41+0.54 

The vast majority of low-mass exoplanet detections to date are 
"hot", involving planets within ~ 1 AU of their host star CuiTently 
only 8 "cool" exoplanets have been detected with masses below 
30 and host separations above 1 AU. This reflects the fact that 
such exoplanets are highly demanding targets for both the transit 
and radial velocity detection methods, techniques which dominate 
current exoplanet statistics. 

Mapping the cold exoplanet regime is crucial for testing and 
informing leading theories of planet formation, such as the core 
accretion and disk instability scenarios. In the core accretion sce- 
nario JSafr onov 1969; Mizuno 1980; Lissauer 1987), planets form 
out of a thick disc of gas and dust by the gradual build-up of ma- 
terial from dust grains into larger objects through collisions. Once 
the objects become large enough, they begin to accrete dust and 
gas via gravity. In the core accretion model, terrestrial planets can 
be considered as the cores of planets that fail to reach the mass 
required for runaway gas accretion, either due to their location in 
the disc or the influence of other planets nearby that grow more 
rapidly. The core accretion process is most efficient in a region 
of enhanced disc density wh ere water and other hydrogen com- 
pound s condense to form ice ( lHavashilll98ll ; IStevenson & Luniri3 
Il988l) . This region (the so-called ice- or snow-line) lies at orbital 
radii ~2.7 AU, likely with only a weak host-mass dependence, 
and is thought to be where most planets form . In the dis c in- 
stabihty scenario ( lKuipei|[l95TI ; ICameronlll978l ; lBosslll997h , gi- 
ant planets form through a gravitational instability in a gaseous 
disc. Disc instability may be th e only mechanism by which gi- 
ant planets can form ( IBoss||20T l|), whilst terrestrial planets ar e still 
thoug ht to form through a process similar to core accretion (IBossI 
l2006h . Migration of planets during formation, due to interactions 
with the disc, can cause both inward taoldreich & Trema inelll980l ; 
IWarj|l997h and outward migration (Masset & Snellgrovel l200lh . 
More violent planet-planet interactions may result in planets being 
scattered inwards (Nagasawa et al. 20 08), outwards or even being 
ejected completely from their systems dVeras et al.l2009b . Tentative 
evidence of unbound (free-floating) planetary-mass objects suggest 
that mor e than one Jupiter -mass planet per star may be ejected in 
this way jSumi et alj201lh . 

Of the relatively few cool low-mass exoplanets detected to 
date at host separations above 1 AU and mass below 30 Earth 
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masses, h alf have been found using the gravitational mi crolensing 
technique jMao & Paczvhskill99l| - lGould & Loebll992h . The peak 
sensitivity of microlensing occurs around the location of the snow 
line, making it a particularly powerful probe of planet formation. It 
is also sensitive to free-floating planets whose existence may pro- 
vide an additional 'smoking-gun' signature of the planet formation 
process. 

Whilst all microlensing surveys to date have been ground 
based, a survey conducted from space is needed to truly open 
up the cold exoplanet parameter space. The probability of a de- 
tectable planetary signal and its duration both scale as proportional 
to Mp, but given the optimum alignment planetary signals from 
low-mass planets are still quite strong. The lower mass limit for 
planets to be detectable via microlensing is reached when the plan- 
etary Einstein r adius becomes smaller than the projected radius of 
the sourc e star (iBennett & Rhielll996h . The ~5.5-A/0 planet de- 
tected by teeaulieu et al.ll2006l) is near this limit for a giant source 
star, but most microlensing events have G or K-dwarf source stars 
with radii that are at least 10 times smaller than this. In order to 
extend the sensitivity to Earth mass and below, it is critical to be 
able to monitor these small source stars that are unresolved from 
the ground. The ideal machine is a wide field imager in space with 
sub arsecond imaging capability. 

The advantages of undertaking a microlensing exoplanet sur- 
vey from space (also discussed in Section [3) were first high- 
lighted some time ago by the study of lBennett & Rhiel ( l2002h who 
looked at the potential science from the Survey for Terrestrial Exo- 
planets (STEP) and Galactic Exoplanet Surx'ey Telescope (GEST) 
mission proposals. Building on these proposals, the Microlensing 
Plane t Finder was propose d to the NASAs Discovery Program in 
2006 teermett et alfcOlOa) . Having realized the synergies between 
the requirements for cosmic shear measurement and microlensing 
planet hunting, a microlensing program was proposed as an ad- 
ditional survey as part of the Legacy Science of the Dark UNi- 
verse Explorer (DUNE) submitted to ES A Cosmic Vision in 2007 
( lRefregieil2 009';'Refr egier et al.l2010l) . Ever since, dark energy and 
microle nsing have bee n advo cated for in a joint mission with white 
papers jBeaulieu et aLll2008|) and at internat i onal c onferences and 
within the community iB eaulieu et ai]|20ld.[20Tll) . Our objective 
is to do a full statistical census of exoplanets down to the mass of 
Mars from free floating to the habitable zone in complements to the 
census from the Kepler mission. DUNE got rebranded into Euclid 
and has been selected as the ESA M2 mission in October 2011, 
with a statistical census on exoplanets via microlensing being part 
of the proposed additional survey in the legacy science. 

The idea promoted in Europe since 2006 of using a single 
space telescope to conduct both a weak-lensing dark energy sur- 
vey and an exoplanet microlensing survey was also f ollowed up in 
the US in a number of white papers and conferences (i Bennett et alj 
I2OO9, 2010a; Gaudi et al. 2009). In US the Exoplanet Task Force 
report (Lunine et al. 2008) to the Astronomy and Astrophysics Ad- 
visory committee concluded that "Space-based microlensing is the 
optimal approach to providing a true statistical census of planetary 
systems in the Galaxy, over a range of likely semi-major axes". 
Following this the US Astronomy 2010 Decadal Survey endor sed a 
comb i ned approach when it top ranked WF IRST (Blandfo rd et alJ 
|20iqj; iBarrv et all 1201 ll ; iGreen et al.l 1201 iL 12012; Dressier et alj 
2012). The subsequent report on impleme nting recommendations 
of the Decadal Review fep ergel et alj|2012[) acknowledges that Eu- 
clid is also capable of undertaking an exoplanet microlensing sur- 
vey. 

Whilst dark energy studies represent the core science of Euclid 



© 0000 RAS, MNRAS 000, 000-000 



ExELS: an exoplanets survey with Euclid I. 



3 



it also aims to undertake other legacy science. The possibility of an 
exoplanet survey is mentioned explicitly in the Euclid Red Book 
jLaureiis et aljEoTlh and is currently under study by the Euclid 
Exoplanets Working Group. This paper presents a baseline design 
for the Exoplanet Euclid Legacy Survey (ExELS). The design for 
ExELS is being developed using a detailed microlensing simulator, 
MAB/xLS, which is also presented in this paper. We focus our atten- 
tion in the present study exclusively to how ExELS will probe the 
cold exoplanet population through microlensing, but ExELS will 
also be able to detect hot exoplanets and sub-stellar objects through 
their transit signatures. This h ot exoplanet science is explored sep- 
arately in a companion paper jMcDonald et alj|20ll hereafter re- 
ferred to as Paper II). ExELS will be the first exoplanet survey 
designed to probe exoplanets over all host separations, including 
planets no longer bound to their host. ExELS will provide an un- 
paralleled homogeneous dataset to study exoplanet demographics 
and to inform planet formation theories. 

We begin the paper by outlining the conservative approach we 
take to our estimates. In Section [3] of this study we overview the 
basic theory behind exoplanet detection with microlensing and we 
also describe the Euclid mission and its primary science objectives. 
In Section|4]we introduce our microlensing simulator (MAB/xLS), 
we describe the Besanjon population synthesis model Galaxy used 
to generate artificial microlensing events, and we also outline a 
baseline design for ExELS. Section|5]presents the results of a sim- 
ulation of the baseline design for ExELS and Section |6] considers 
the effects of a number of variations and de-scope options to the 
baseline design. We end with the summary discussion in Section|7] 



2 A CONSERVATIVE APPROACH 

Throughout our study of the capability of Euclid for detecting ex- 
oplanets we adopt a conservative approach. There are two reasons 
for this. Firstly, the design of Euclid is itself still evolving. Sec- 
ondly, since time on a space telescope is expensive, a feasibility 
study such as carried out in this paper must demonstrate that key 
science goals are likely to be achieved rather than merely being an 
aspiration. 

This means that, wherever possible, we aim to make detailed 
predictions anchored to models which are known to agree with cur- 
rent data. Where details of models require some assumptions these 
assumptions must not be overly-optimistic. An example of this ap- 
proach is our simulation of photometry. The most accurate repre- 
sentation of the photometric methods that will be used on the real 
data would be to simulate PSF fitting or weighted aperture pho- 
tometry. However, crowded field photometry is notoriously diffi- 
cult, and there will always be cases where automated data analysis 
pipelines will fail to perform the photometry optimally. Simulating 
all the possible complications in the photometry is impossible. If 
we were to simulate PSF fitting or weighted aperture photometry, 
complications that degrade or destroy the photometry would not be 
modelled and the assumption we had made, while being accurate 
would be optimistic. Instead, we choose to simulate the photometry 
as unweighted aperture photometry (see section l431 for full details). 
Aperture photometry is a less accurate representation of the actual 
data analysis methods that will be used, and is less effective than 
optimal photometry by a small but sometimes significant amount. 
However, this choice is conservative, and helps to ensure that we 
do not overpredict the performance of the mission. 

Given our conservative approach we can have confidence that 
the scientific yields we predict are realisable with Euclid. 



3 EXOPLANETARY MICROLENSING FROM SPACE 

Gravitational microlensing describes the transient deflection and 
distortion of starlight on milli-arcsecond scales by interv ening 
stars, stellar remnants or planets (for a recent review see IMad 
Microlensing is distinguished from ordinary gravitational 
lensing in that whilst multiple images are produced, they are not 
resolvable. Instead one observes a single apparent source which 
appears magnified by a factor 

Uy/u^ + 4 

where the impact parameter m is a dimensionless angular separa- 
tion between the lens and source measured in units of the angular 
Einstein radius of the lens. A microlensing event is observable as a 
transient achromatic brightening of a background source star last- 
ing for ~ 6 — 60 days, where tE is the Einstein radius cross- 
ing time. For a single lens the lightcurve profile is time- symmetric, 
with a peak magnification occurring when the impact parameter is 
at its minimum u = uq. The lensing signal from foreground stars 
is detectable in a few out of every million background stars located 
in crowded stellar fields such as the Galactic bulge. A planet orbit- 
ing a foreground lensing star may, in a few percent of microlensing 
events, perturb the microlensing signal causing a brief deviation 
which lasts for tp ~ Afp/M,, where M« and Mp are the host 
and planet masses. Typically tp is in the region of a day for a Jupiter 
mass planet down to a few hours for Earth mass planets. The in- 
trinsically very low probability ~ 0(10^*) of an exoplanetary mi- 
crolensing signature against a random background source star, cou- 
pled with the brief deviation timescale associated with Earth-mass 
planets, places huge technical demands on microlensing surveys. 

The probability of a planetaiy perturbation occurring scales 
roughly as the square root of the planet mass, or more strictly, as th e 
square root of the planet-host mass ratio q jGouId & Loeb|[l992h . 
This shallow sensitivity curve makes microlensing ideal for detect- 
ing low-mass planets. The scaling breaks down below about a Mars 
mass, where finite-source effects begin to w ash-out planetary signa - 
tures, even for main-sequence source stars dBennett & Rhiell2002l) . 

The sensitivity of microlensing to planets peaks close to the 
Einstein radius te with projected semimajor axis ax ~ te ~ 
2 AU, corresponding to where the microlens i ng images are most 
likely to be perturbed (Wa mbsganssI 1 19971 . iGriest & Safizadehl 
Il998h . However there is significant sensitivity to planetary orbits 
wi th Ox ~ 0.5 AU , and outwards to i nfinity (i.e. free-floating pldin- 
ets lHan et alJl2004l ; ISumi et alj|201 ih . 

Owing to its high stellar density and microlensing optical 
depth, the Galactic bulge is the best target for microlensing stud- 
ies. Towards the bulge, extinction is a significant problem at op- 
tical wavelengths. Additionally, the extreme stellar crowding and 
arcsecond-scale seeing mean that only the giant star population 
can be properly resolved from the ground (iBermetd 12004)) . Ob- 
serving in the near-infrared lessens the eff ects of dust and so pro- 
vides a larger microlensing optical depth dKerins et al.ll2009l) but, 
from the ground, stellar crowding problems are even more severe, 
and noise levels are enhanced due both to the sky and unresolved 
stellar backgrounds. Therefore, in order to monitor enough source 
stars, ground-based surveys must regularly cover ~100 deg^. 
Curre nt and futur e ground-based surveys - e.g.. MOA-II dSumil 
|20IO|) . OGLE-IV (Udalski 20TT), KMTNet jKim et alll20I ()) and 
AST3 l lYuan et a l. 2010) - with wide-field imagers will achieve 
suitable cadence and areal coverage to detect routinely large num- 
bers of giant planets if they exist in sufficient abundance. How- 
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ever they will not be able to monitor enough stars at high-cadence 
to detect Earth-mass planets at a significant rate. For this reason, 
targeted follow-up of promising microlensing events by large net- 
works of small telescopes is currently used to achieve high ca dence 
and continuous event coverage (see, e.g.. lGoulderZll2010h . and 
to push the sensitivity of ground-based microlensing firmly int o 
the super-Earth regime (iBeaulieu et al ] |2006t iBennett et al.ll2008h . 
However, the follow-up networks only have the capacity to ob- 
serve ~100 events per year or fewer with sufficient cadence (IPealel 
l2003h . This allows the mass function to be probed down to ~5- 
1OM0, and possibly the semi-major axis distribution of planets 
above ~5OM0, but is unlikel y to provide more than isolated de- 
tections below these masses l lPealell2003l : lBennetJ|2004l : iDominikI 

(mil). 

Observations from space are able to overcome many of the 
problems facing ground-based observers. A space telescope has 
better resolution due to the lack of atmosphere and also a lower 
sky background, especially in the infrared. This means that with 
appropriate instrumentation, a space telescope can resolve main- 
sequence sources in the bulge and monitor the required ~10^ 
sources over a much smaller area. This in tu rn allows high-cadence 
observations on a small number of fields jBennett & Rhiell2002l : 
lBermett|[2004l) . The fundamental requirements of a space telescope 
for a microlensing survey are a wide field of view (> 0.5 deg^), 
with a small pixel scale. In order to minimize the effect of extinc- 
tion towards the Galactic bulge, it should observe in the near in- 
frared. The telescope must also have a large enough collecting area 
to allow high-precision photometry of main-sequence bulge stars 
in short exposure times. These are almost exactly the same require- 
ments as for dark energy studies using weak lensing, which are 
already driving the hardware design of Euclid. 

3.1 The Euclid mission 

Euclid is an M-class mission within the ESA Cosmic Vision pro- 
gramme. It aims to investigate the nature of dark energy through 
measuremen ts of weak gravitatio nal lensing and baryon acoustic 
oscillations jLaureiis et alj l201~ll) . Euclid will comprise a 1.2-m 
Korsch telescope with a high-resolution optical imager (VIS) and 
a near infrared imaging spectrometer (NISP), operating simultane- 
ously. The core science mission will involve a 15000-deg^ wide 
survey and 40-deg'^ deep survey over six years to measure galaxy 
shapes and photometric and spectroscopic redshifts. VIS will ob- 
serve with a wide optical band-pass covering R, I and Z, and NISP 
will have available three infrared filters: Y , J and H . The currently 
envisaged step and stare survey strategy of Euclid means that for 
up to two months per year it will point towards the Galactic plane 
and aw ay from its primary science fields. As stated in lLaureiis et alj 

l|) it is intended that some of this time will be devoted to other 
legacy sci ence. A planetary mi crolensing survey is one option de- 
scribed in lLaureiis et al.l ( 1201 ih and is being actively evaluated by 
the Euclid Exoplanets Working Group. 

The similarity of hardware requirements for dark energy and 
exoplanet microlensing space missions has been recognised for 
some time (Bennett & Rhie 2002), and most recently by the 2010 
US Astrophysics Decadal Review dBlandford et al.l |20T(]|) . This 
review recommended the merger of three mission concepts into 
one mission, the W ide-Field Infrared Survey Telescope ( WEIRST, 
lOreen et"H]|2012h . Two of the core science objectives for WEIRST 
are a dark energy survey and an exoplanets survey using microlens- 
ing. In the baseline WEIRST concept the microlensing survey will 
total 432 days, somewhat longer than will be feasible for ExELS. 



4 THE MANCHESTER-BESANCON MICROLENSING 
SIMULATOR (MAB^iLS) 

We have designed the Manchester-Besanfon micro-Lensing Sim- 
ulator (MAB/^LS - pronounced raaj-buls) to perform detailed sim- 
ulations of the ExELS concept. MAB/iLS is the first microlensing 
simulator to use a combination of a population synthesis Galactic 
model with a realistic treatment of imaging photometry. This means 
that every aspect of the simulation, including the event rate calcu- 
lations, blending and photometry are simulated self-consistently. 

Several key ingredients are needed in order to simulate any 
microlensing survey. A simulator must draw its simulated events 
from a Galactic model and distributions of the event parameters. It 
must simulate the observations of the survey, and finally, it must 
also simulate the detection criteria used to select its sample of 
events. It is also necessary to make a choice as to the complexity 
of the microlensing model used to simulate events. For example, is 
the lens composed of a single mass or multiple components? Are 
higher-order effects such as parallax and orbital motion included? 
In the rest of this section we will discuss both how MAB/iLS imple- 
ments each component of the simulation and the choice of param- 
eters we use in the simulation of ExELS. Unless stated otherwise, 
we have taken the s urvey parameters from the Euclid Red Book 
dUaureiis et al.ll2oTll) . 

4.1 The Besan^on Galactic model 

Underpinning the MAB/.tLS mic rolensing event generation is th e 
Besangon model dRobin & Crez e 1986; Robin et al. 2003.l 2012h . 
a population synthesis model of the Galaxy. The Besangon model 
comprises five main stellar populations, a spheroid (stellar halo), 
thin and thick discs, a bar and bulge. The stars of each popula- 
tion are assumed to be formed from gas for input models of star 
formation history and initial mass function (IMF). The stars are 
aged according to model ev olutionary tracks to their present-day 
state iHavwoodetalJIiggTl) . This determines the distribution of 
stellar bolometric fluxes, which are converted to colours and mag- 
nitudes using stellar atmosphere models convolved with standard 
band-pass templates in various photometric systems. 

The spatio-kinematic distribution of the disc stars is deter- 
mined by integration of a self-consistent gravitational model us- 
ing the Poisson and Boltzmann equations. Finally, the observed 
colours and magnitudes ar e corrected for extinction using a three- 
dimensional dust model cMarshall et al.ll2006) . A limited number 
of model parameters are then optimized to reproduce observed star 
counts and kinematics. The output of the model is an artificial cata- 
logue of stellar photometry and kinematics for a survey of specified 
sensitivity and areal coverage. 

The Besango n model is in constant development (e.g., 
iRobin et alj |20I2|) . In this work we use version 1106 of the 
Besangon model, though an updated version of the model has been 
released since. In subsequent models, the properties of the Galac- 
tic bar and bulge (see below) change somewhat from those we use 
here. Below we briefly overview the properties of the main stellar 
components used to generate microlensing events in MAB/iLS. The 
Solar Galacto-centric distance in the model is 8 kpc. 

4.1.1 The .stellar halo 

The stellar halo is modelled as being formed by a single burst of star 
formation at 14 Gyr, with metallicity centred at [Fe/H] — —1.78 
and with a dispersion of 0.5. It has a triaxial velocity distribution 
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with dispersions (o-y, crv, ow) ~ (131, 106, 85) km s^^. Its den- 
sity is small everywhere, even at the Galactic center, and so con- 
tributes only marginally to the optical depth and microlensing event 
rate. 



4.1.2 The bar 

The bar, a l tered from the bulge-like component used by 
iKerins et alj j2009t) . co nsists of a boxy t riaxial distribution, simi- 
lar to that described by IPicaud & RobinI ( 2004 ), bu t with a Gaus- 
sian density law as o pposed to a FreudenreichI ( Il998[) sech^ 
law ( lRobinetalJl2012h . The major axis of the triaxial structure 
lies at an angle of 12.5° relative to the Sun-Galactic centre line 
of sight and has scale lengths (X, Y, Z) = (1.63, 0.51, 0.39) kpc, 
where the X direction is parallel to the major axis and the X and 
Y axes lie in the Galactic plane. It is truncated at a Galactocentric 
radius of 2.67 kpc. The bar rotates as a solid body with a speed 
of 40 km s~^ kpc~^. The velocity dispersions in the bar along the 
axes defined above are (113, 115, 100) km s"\ The central stel- 
lar mass density of the bar, excluding the central black hole and 
clusters, is 19.6 x lO'-'M© kpc"^. 

Embedded within the bar is also an additional compo- 
nent (somewhat d ifferent from the "thick bulge" component in 
iRobin et al.|[2012l) . However, in the version of the model we use 
here, its density is smaller by ^10"* times that of the bar, so we 
do not describe it further. We use the terms "bar" and "bulge" inter- 
changeably from here onwards to mean the bar component of the 
Besanjon model. 



Table 1. The magnitude range, and the average number of stars (A'^*) in the 
Besan^on model star catalogues used in this work. Bright, moderate and 
faint are the three levels of catalogues used to build the simulated images. 



The stellar po pulation of th e bar is assumed to fo rm in a single 
burst 7 .9 Gyr ago jPicaud & R obin 2004), following TGirardi et all 
( I2OO2I) . The bar IMF (diV/dAf) scales as A/"^ between 0.15Mq 
and O.7M0 , and follows a Salpeter slope above this mass. The pop- 
ulation has a mean metallicity [Fe/H]= 0.0 with dispersion 0.2 and 
no metallicity gradient. The stellar luminosities are calculated us- 
ing Padova isochrones CGirardi et al..200Z) . 



4.1.3 The th ick disc 

The thick disc is modelled by a single burst of star formation at 
11 Gyr. Its properties have been constrained using star counts by 
iRevle & Roblrj j200ll) . The thick disc contributes only marginally 
to the microlensing event ra te, so we do n ot describe it in detail. Its 
parameters are described bv lRobin et all [2003i) . 



4.1.4 The thin disc 

The thin disc is assumed to have an age of 10 Gyr, over which 
star formation occurs at a constant rate. Stars are formed with a 
two-slope IMF that scales as a power-law M~^'^ from O.O79il/0 
to IMq and above, based on the Hipparcos luminos- 

ity fiincti£n_(e.gjHOTWood^^ with updates described 

by iRobin et"aL j2003l). Stars be l ow 1M b follow the evolution- 
ary tracksofv^ndenBerg'^r^ ( I2OO6I) , while those ab ove fol- 
low Schaller et al. ( 1992) tracks. The thin disc follows an lEinastol 
( IT979 ) density profile with a central hole. The density normal- 
ization, kinematics and metallicity distribution of the disc de- 
pend on stellar age, with seven a g e ran ges defined, whose pa- 
rameters are given by iRobin et"ar I ( I2OO3I) . The Solar velocity is 
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the youngest disc component which has disc and hole scale lengths 
of 5 kpc and 3 kpc, respectively. The disc is truncated at 14.0 kpc. 
The scale height of the disc is computed self-consistently using 
the Galactic potential via the Boltzmann equation as described by 
iBienavme et al. jl987l) . Also modelled in the disc are its warp and 
flare ( iRevle et alj|2009l) . 



4.1.5 Extinction 

Extinction is computed using a three-dimensional dust distribution 
model of the i nner G alaxy < 100°, |6| < 10° ), built by 
'Mars hall et"ZI (I2OO6) from analys is of 2MASS data (ICutri et alj 
,20031) using the Besanjon model. iMarshall et"al] ( |2006|) did this 
by comparing observed, reddened stars to unreddened simulated 
stars drawn from the Besangon model. From this the extinction as 
a function of distance along a given line of sight is computed by 
minimizing between observed and simulated J — Ks colour 
distributions. The resulting map has a ~15-arcmin resolution in I 
and fe, and a distance resolution ~0. 1-0.5 kpc, resulting from a 
compromise between angular and distance resolution. 



(f/0, V0,Wq) = (10.3,6.3,5.9) km s-\ with respect to the 
local standard of rest Vlsr = 226 km s-^. The disc has a scale 
length 2.36 kpc, and the hole has a scale length 1.31 kpc, except for 



4.1.6 Other components 

The Besanjon model also takes account of other Galactic com- 
ponents, including the mass due to the dark matter halo and in- 
terstellar medium. The details of these components are given by 
'Robi n et alj j2003l) . White dwarfs are included in the model sepa- 
rately to normal stars, with separate densities and luminosity Junc- 
tions determined from observational constraints (IRobin et alj2003l 
and refere nces therein). The evo lutio nary tracks and atmosphere 
models of [Bergeron et al.l d 1 995b and lChabrieil ( Il999l) are used to 
compute their colours and magnitudes. 



4.2 Microlensing with the Besan^on model 

Following the method of iKerins et alj ^20091) , MAB/xLS uses two 
star lists output by the Besangon simulation to construct catalogues 
of possible microlensing events and calculate their properties. The 
first list, the source list, is drawn from the Besangon model using a 
single magnitude cut in the primary observing band of the survey. 
A second list, the lens list, is drawn from the model without a mag- 
nitude cut. Both source and lens lists are truncated at a distance of 
15 kpc to improve the statistics of nearer lenses and sources that 
are much more likely to be lensed/lensing. 

Overall microlensing event rates are calculated along mul- 
tiple lines of sight, with spacing set by the resolution of the 
[Marshall et Zl ( I2OO6I) dust map. The total rate due to each pair of 
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source and lens lists, about the line-of-sight [l, b), is 



„ Sources / ^ Lcnsc 

Silos I 



(2) 



where flios is the solid angle covered by a dust-map resolution- 
element, and S^ls and SQi are the solid angles over which the 
source and lens catalogues are selected, respectively. The rate is 
calculated over all possible source-lens pairs to minimize the noise 
of counting statistics. The average sizes of all the catalogues used 
in this work are listed in Table [T] 

To simulate microlensing, MAB/iLS draws sources and lenses 
from their respective lists with replacement, requiring the source 
be more distant than the lens. From the source and lens parameters, 
the Einstein radius and timescale are computed, as well as the rate 
weighting assigned to the event 



7 = ""On 



(3) 



where Momax is the maximum impact parameter of the event; how 
womax is determined is discussed in the following sections. Events 
are simulated and those that pass the detection criteria are flagged. 
The rate of detections in a given dust-map element is the sum of 
the weights of detected events normalized to the sum of the rate 
weightings for all the simulated events ~ this is essentially a de- 
tection efficiency. The detection efficiency is then multiplied by the 
total line-of-sight rate computed in Equation|2]to yield the expected 
detection rate for 0.25 x 0.25 deg^, the size of the dust-map ele- 
ment. These rates are then summed over all the dust-map elements 
to yield the total simulation event rate. 

For the ExELS simulations we restrict the source magnitude to 
a vega-based iif-band magnitude J/vcga < 24. This corresponds to 
an AB magnitude limit of Hab < 25.37. Unless otherwise noted, 
AB magnitudes will be used throughout the paper. 



4.2.1 Normalization of the event rate 

The absolute number of stars in the Besan^on model has been set 
by normalizing their number density to match star counts along a 
number of l ines of sight that sample the different components of 
the Galaxy jRobin et alj|2003h . This process will only be accurate 
down to the limiting magnitude of the star count data, and it is 
possible that the number of fainter stars predicted by the model 
could be incorrect. While these fainter stars may not contribute to 
star counts, they do contribute to the microlensing event rate, either 
as lenses or sources. It is therefore important to make sure that the 
microlensing event rates predicted by the Besanfon model match 
those that are observed. 

The microlensing event rate is a potentially powerful but com- 
plex probe of Galactic structure because it depends on the kine- 
matics of the lens and source populations as well as the lens mass 
function. Often, surveys aim instead to measure the microlensing 
optical depth, which is much more cleanly defined as it only de- 
pends only on the distribution of lenses and sources along the line 
of sig ht. T he total event rate T is related to the optical depth r by 
(e.g., lPaczvhski.,1996^ 



r oc 



(4) 



where A?* is the number of monitored sources and {I-e) is the av- 
erage event timescale. We can verify the microlensing event rate 
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Fig ure 1. Comparison of the luminosity function in the bulge as measured 
by Hol tzman etaP il998h in Baade's window {t=l,b= -3.9) to that of 
the Besangon model at the same location. The measured luminosity func- 
tion, shown by the black line, has been returned to the /-band apparent 
magnitude scale by adopting t he distance mod ulus {fi = 14.52) and ex- 
tinction {Aj = 0.76) values of lHoltzman et al.l fl998i. The red line shows 
the Besangon model's luminosity function at Baade's window, while the 
blue line shows the Besangon model luminosity function at the center of 
the ExELS fields (£ = 1.1. = —1.7). The grey line shows the Earth- 
mass planet detection efficiency for Euclid as a function of source /vega 
magnitude (abitrarily shifted on the log scale). 



predicted by the Besanfon model by comparing each of the quan- 
tities on the right hand side of Equation |4]to measured values, and 
make a correction to the rate if required. 

In the bulge region, the Besanjon model has been normal- 
ized to star coun ts from the 2MASS survey dCutriet al.lEool 
iRobin etalj|2012h . which is relatively shallow compared to the 
sources that Euclid will observe. The number counts of fainter 
sources are extrapolated using the IMF and extinction maps, and 
any uncertainties in these will propogate to the source counts. 
There is a relative paucity of deep star count measurements in the 
bulge with which to test the Besanfon model assumptions, with 
published measurements only along a single line of sight close to 
our proposed Euclid fields. To assess the validity of TV*, we com- 
pare t his measurement of th e luminosity function in Baade's win- 
dow jHoltzman et alj[l998l) to the Besangon star catalogue pro- 
duced for the same line of sight. Figure [T| shows both luminosity 
functions. There is good agreement between the two luminosity 
functions in the magnitude range 17 < /vcga < 20, but fainter 
than this the Besangon model under- pr edicts the numb e r of s tars. 
Integrated over the whole range of the iHoltzman et al. IdHii lu- 
minosity function, the Besangon model predicts 32.46 x 10^ stars 
per square degree, but Holtzman et al. ( 1998) measure 42.06 x 10*' 
stars per square degree. To correct the event rate for this lack of 
stars requires multiplying by a factor of 1.30. While we adopt this 
correction, we caution that a significant fraction of the discrep- 
ancy arises from the faintest part of the luminosity function, where 
Holtzman et al. ( 1998) argue that their co mpleteness corrections 
are uncertain. Beyond the faintest bin of the lHoltzman et all ( 1 1 998b 
luminosity function we might worry that the number of stars keeps 
on increasing in reality, while in the Besangon model it begins to 
fall off, suggesting a larger correction factor would be required. 
However, average extinction in the ExE LS fields (Aj — 1.73 at 
a distance of 8 kpc [Marshall et al.ll2006l) i s nearly one magnitude 
more than Ai = 0.76, the value adopted bv lHoltzman etaljjl998h 
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Figure 2. A comparison of measured optical depths to red clump giants to 
those calculated from the Besangon model. Open square, filled circle and 
asterix data points show results from the MACHO (Popowski et al. 20(3), 
OGLE <Sumi et al. 2006) and EROS iHamadache et al. 2006) surveys, re- 
spectively. The solid hne shows the average optical depth to red clump stars 
selected from the Besanfon model based on their absolute magnitudes and 
colours. The dashed line shows the same line as the solid curve, but multi- 
plied by a constant 1.8 to match the data. 



in t heir Baade's window fie ld. This implies that in the ExELS fields 
the iHoltzman et alj ( Il998 ) luminosity function extends down to 
/vega ~ 25.2. At this source magnitude, the planet detection ef- 
ficiency for Euclid has fallen to roughly a third of the average at 
brighter magnitudes, and falls rapidly as the sources get fainter. 
Therefore, while there may be more faint stars that the Besanjon 
model is missing, including these source stars will not significantly 
increase the number of planet detections. 

The optical depth in the Galactic bulge has been measured 
to two different source populations: red clump giants and dif- 
ference imaging analysis (DIA) sources. Measurements of DIA 
optical depths are typically slightly higher than those for clump 
giant sources. Clump giants are abundant, bright standard can- 
dles, and those in the bulge can be easily recognised and isolated 
by their position on a colour-magnitude diagram. They therefore 
make an ideal tracer population of the bulge. DIA sources how- 
ever, are less clearly defined, as they depend on the sensitivity 
of the survey. Due to this difficulty of defining the DIA source 
population, we only compare the Besanjon model to measured 
red cl ump optic al depths, which hav e been measured b y the MA- 
CHO <PoDOWski et al. 2005) . EROS faamadache et aljf2006) and 
OGLE jSumi et al]|2006l) survevs. Figure |2] shows the red clump 
optical depth measurements of each of these surveys as a func- 
tion of absolute Galactic latitude, together with the average opti- 
cal depth to red clump stars as predicted by the Besangon model. 
The Besangon optical depths are averaged over the longitude range 
—0.525 < I < 2.725. It is clear that the Besanjon optical depth is 
somewhat lower than the measured optical depth. Multiplying the 
Besangon model optical depths by a constant factor of 1.8, brings 
the predictions into good agreement with the measurements. 

The average timescale reported by microlensing experiments 
is somewhat ill-defined due to the arbitrary timescale cut-offs ap- 
plied to their event samples, but are typically in the range of 20- 
30 days. The average timescale calculated from the timescale dis- 
tribution presented in Figure[T3]is 21.4 days without applying any 
detection criteria and 29.2 days for events detected above base- 



line. The average timescale of the ISumi et al] ( 1201 ih sample is 
26.0 days . The sample det ected above baseline is most compara- 
ble to the dSumi et al I20T l|) sample, but is not directly comparable. 
As the difference is of the order of 10-20 percent, and the samples 
are not directly comparable, we choose not to make a correction to 
the event rate based on the timescales. 

Combining the required correction factors for optical depth 
and source counts, we conclude that the microlensing event rates 
predicted by the Besangon model likely require a correction factor 



/1106 = 1.8 X 1.30 = 2.33, 



(5) 



where the subscript emphasizes the fact that this scaling is only 
applicable to the version of the Besangon model we are using. We 
have multiplied all raw results by this factor throughout the paper. 
To account for further uncertainty in the overall event rate we define 
an additional multiplicative factor /r, with which all the event rates 
and numbers of planet detections should be multiplied. In this paper 
we advocate for a value /r = 1. Note, however, that dCreen et al.l 
[2012) choose a value /r = 1.475 to compromise between different 
values of the planet detect ion efficiency de t ermin ed from MABELS 
and simulations following jBennett & Rhid j2002l) . 

One possible cause of the low predicted optical depth could 
be missing low-mass stars and sub-stellar objects too faint to 
be included in star counts. The lower cut-off of the bulge 
mass function in the Besangon model is O.15M0. Extending the 
mass function down to ~O.O3M0, keeping the same low-mass 
slope (M~^), would acc ount for the missing optical depth (see 
ICalchi Novati et aT 2008, for a discussion of the effect of the mass 
function of microlensing event rates). Adding such low-mass stars 
to the star catalogues would increase the number of planet detec- 
tions by a factor larger than the increase in optical depth, because 
the mass ratio of planets around these stars would be larger than 
the mass ratio of planets around a star of the average stellar mass 
in the catalogues. Another possible cause of the low optical depth 
is the lack of high-mass stellar remnants - neutron stars and black 
holes, which are not included in the Besangon catalogues. Should 
these be the cause of the low optical depth, the number of planet 
detections would not increase as much as the optical depth, be- 
cause even if planets remained around these objects, their mass 
ratios would be smaller than the average in the unmodified cata- 
logues. Other possible causes of the optical depth discrepancy, such 
as problems with Galactic structure, would cause the number of 
planet detections to change in the same way as the microlensing 
event rate. Note that in all the scenarios discussed here, the av- 
erage timescale is affected - adding low-mass stars decreases the 
average timescale, giving a further boost to the event rate, while 
adding high-mass remnants increases the average timescale, further 
supressing any boost. Comparing the Bes angon timescale distribu- 
tion to that observed bv lSumietaLl ( l201lh (see Figure [13} suggests 
that the Besangon model is missing short timescale events, but each 
timescale distribution has different selection criteria, so it is impos- 
sible to draw firm conclusions. 



4.3 The microlensing events 

MABELS uses user-supplied functions to compute microlensing 
lightcurves including any effects that the user wants to model. 
For this work we modelled only planetary lens systems composed 
of a single planet orbiting a single host star. As we want to in- 
vestigate the planet detection capability of ExELS as a function 
of planet mass Mp and semimajor axis a, we chose to simulate 
systems with various fixed values of planetary mass in the range 
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0.03 - 10" M0 and semimajor axis distributed logarithmically in 
the range 0.03 < a < 30 AU. We assume a circular planetary orbit 
that is inclined randomly to the line of sight. The orbital phase at 
the time of the event is again random; we do not model the effects 
of orbital motion in the lens. The impact parameter and angle of 
the source trajectory are distributed randomly, with the impact pa- 
rameter in the range uo = 0-itomax. For ExELS simulations we 
choose uomax = 3. While it may be the case that some of the stel- 
lar microlensing events with « 3 will not be detected, it is still 
possible for planets to cause detectable signals. 

The planetary microlensing lightcurves are computed assum- 
ing that the source has a uniform intensity profile (in other words, 
no limb darkening). Test simulations including the effect of limb 
darkening (which is small in the infrared) showed that inclusion 
of the effect would change the number of Earth-mass planet de- 
tections by less than 1 percent. The finite-source magnification 
is computed using the h exadecapole approximation when finite - 
source effects are small (Peicha & Hevrovskvl [2009'; 'Gould"2008^ 
and t h e contou ring method when they are not ([Gould & Gaucherel 
ll997l ; fDominiklll998 ). Finite-source effects are accounted for in 
single -lens lightcurve calculations using the method o flWitt & Mad 
( I1994I) . When fitting lightcurves with the single-lens model, we use 
a finite-source single-lens model if the impact parameter uq < lOp, 
where p is the ratio of angular source radius to the angular Einstein 
radius. Otherwise the point-source single-lens model is used. 

4.3.1 Free-floating planets 

Observations from nearby star clusters, as well as tentative evi- 
dence from ground-based microlensing surveys, suggests that plan- 
ets can occur unbound from any host, sometimes referred to as free- 
floating planets. If free-floating planets exist in significant numbers 
then ExELS should detect them as relatively brief single-lens mi- 
crolensing events. 

At this stage we have no clear information to allow us to char- 
acterise a Galactic population of free-floating planets with confi- 
dence. What we know from young star clusters like sigma Orionis 
is that brown dwarfs (0.072 — 0.013 Mq) and massive free-floating 
planets (0.013 — 0.004 Mq) are as numero us a low-mass stars 
with masses in the interval 0.25 — 0.072 Mq jPena Ramirez et al] 
12012 !). However, given the uncertainties over the characteristics of 
a Galactic-wide distribution of planets we adopt a simple scalable 
assumption of one free-floating planet of mass A/ffp per Galactic 
star. 

As free-floating planets are single, point-mass lenses we treat 
them in a separate MAB/iLS simulation. Each lens star drawn from 
the Besanjon simulation is replaced by a planet of mass A/gp. We 
simulate a range of values for Algp from 0.03 — 10* M(^. We 
assume the planets retain the same kinematics, but the fundamen- 
tal microlensing properties such as the Einstein radius change to 
reflect the reduced mass. Ejected planets may well have a some- 
what larger velocity dispersion than their original hosts, in which 
case the rate of free-floating planet events increases proportionately 
and the timescale decreases inversely with their velocity. We as- 
sume that free-floating planets emits no detectable light, which is 
a good assumption for typical distances at which planets are de- 
tectable through microlensing. Each lightcurve is calculated using 
the finite-source single-lens model. The impact parameter is cho- 
sen to lie in the range ito — 0-itomax, where we choose itomax ~ 1 
to remain conservative, and we require that the time of peak mag- 
nification lie within an observing season (unlike for the standard 
simulations). 




Figure 3. The approximate location of the thi'ee ExELS field pointings 
(solid line rectangles) assumed for the simulation, with the middle of the 
three fields centred at / = 1.1°, 6 = —1.7°. Each ExELS field has dimen- 
sions of 0.76 X 0.72 deg2. The background image of the Galactic Centre is 
a near-infrared mosaic of images from the VVV survey jSaito et al.ll2012h . 
Background image credit: Mike Read (Wide-field Astronomy Unit, Edin- 
burgh) and the VVV team. 



4.4 Euclid observing strategy 

The ExELS survey must be capable of detecting planets at least 
down to Earth masses, which means we require an observing ca- 
dence of no more than 20 mins between repeat observations of the 
same field. It must also monitor enough source stars over a suffi- 
ciently long observing baseline to ensure a healthy detection rate. 
As shown bv lKerins etZH2009l) the event rate is optimised at near- 
infrared wavelengths, suggesting that the NISP camera should be 
the primary instrument for ExELS. In Section |6T| we show that 
this is indeed the case, despite the significantly worse resolution of 
the NISP instrument relevant to the optical VIS instrument. 

In order to achieve a cadence of no worse than 20 min, ExELS 
will be able to monitor up to 3 target fields of ~0.5 deg^ with a 
total exposure of 270 s per pointing, split into stacks of 3 (Y- and 
J-band) or 5 (iJ-band) exposures with NISP. We assume that there 
is 5 s of dead time between the exposures of a stack. The VIS in- 
strument pointings consist of a single 540-s exposure. We assume a 
baseline slew and settle time of 85 s, though in Section |63] we also 
consider the effect of a substantially longer slew and settle time. We 
assume that any readout, filter wheel rotation and data down-link is 
performed during slewing. Some of these parameters are summa- 
rized in Table |2] 

We simulate a total observing baseline of 300 days for Ex- 
ELS, spread over 5 years with two 30-day seasons per year. This 
strategy is determined by the design of the spacecraft's sun shield. 
This restricts Euclid to observing fields with solar aspect angles 
between 89 and 120 degrees. As the Galactic bulge lies near to the 
ecliptic, Euclid can only observe bulge fields uninterrupted for up 
to 30 days, twice per year. It should be stressed that a lO-month sur- 
vey represents a firm theoretical maximum that could be possibly 
achieved over 5 years. In practice the cosmology primary science 
will likely prohibit much legacy science being undertaken in the 
first few years of the mission, so that a 6-month exoplanet survey 
likely represents a more achievable goal during the 6-year primary 
cosmology mission. It is possible that, if Euclid remains in good 
health beyond 6 years, a full 10 month programme could be com- 
pleted after the cosmology programme is complete. We therefore 
investigate the exoplanet science returns possible for a survey of 
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up to 10-months total time. We assess the impact of shorter total 
baselines in Section[5] 

For the ExELS simulation we use three contiguous Euclid 
pointings aligned parallel to the Milky Way plane, with the central 
field located at Galactic coordinates I = 1.1°, b = -1.7° (J2000), 
as shown in Figure |3] Each Euclid field covers 0.76 x 0.72 deg^ 
on the sky, giving a total ExELS survey area of 1.64 deg^. We con- 
servatively assume most of the observations are taken with NISP in 
only one filter (we show in Section |6] that _ff-band is the best fil- 
ter choice), at a cadence of roughly one observation every 20 mins. 
Conservatively we add colour information from the two other NISP 
filters and the VIS camera at a rate of only one observation every 
12 hours. This conservative assumption guarantees we will not be 
limited by telemetry rate restrictions. However, in Section [6T| we 
consider the benefits of simultaneous NISP and VIS imaging. For 
the hot exoplanet science investigated in Paper II we note that it is 
important to achieve high cadence observations with both VIS and 
NISP instruments. Therefore strong limitations in telemetry could 
impact somewhat upon the hot exoplanet science but is unlikely 
to impact significantly upon the cold exoplanet science discussed 
here. 

The number of planets which can be detected by ExELS will 
be governed by the overall rate of microlensing within the survey 
area, though only a small fraction of these will have detectable 
planetary signatures. The expected overall number of microlensing 
events (with or without planet signatures) that would be detected as 
significantly magnified by ExELS is ~27000 events with uq < 3 
over the course of a 300-day survey. This is in excess of the total 
number of microlensing events discovered by all ground-based mi- 
crolensing surveys since they started operations twenty years ago. 
However some of these events may not be well characterised if 
they peak outside of one of the observing seasons. Placing the re- 
striction that the time of peak magnification, to, must be contained 
within one of the observing seasons lowers the overall number to 
~9800 events (~5700 events with uo < 1) for a 300-day cam- 
paign, or about 1000 events per month. This is an improvement of 
a factor of ~ 65 in d etection rate pe r unit time per unit area over 
the OGLE-IV survey ( IUdalskill20T ll) in its best field, which yields 
~ 15 events per month per deg^. Between now and Euclid's sched- 
uled launch in 2019, the OGLE-IV survey observing the bulge ~8 
months per year can detect a similar number of microlensing events 
to Euclid observing for 10 months total. However, the ExELS sur- 
vey will be much more sensitive to low-mass planets per event. 



4.5 Photometry 

In order to accurately account for the effects of severe stellar crowd- 
ing on photometry of Galactic bulge stars, MAB/iLS produces sim- 
ulated images for each microlensing event it simulates. The source 
and lens star of each microlensing event are injected into the image, 
with the source star's brightness updated at each epoch. Finally, 
relative aperture photometry is performed to measure the source 
brightness in each image. 

The image is constructed from a smooth background compo- 
nent and stars drawn from the Besan^on model catalogues. Stars are 
added to the images at random locations on a fine (9 x 9) sub-pixel 
grid, using numerical PSFs that have been integrated over pixels. 
Each star is tracked so that it is included at the correct position and 
brightness in images taken with different filters or instruments. In 
this way, blending is computed consistently throughout the simula- 
tion. In order to avoid small-number statistics for bright stars with- 



out using huge catalogues, we use tiered catalogues with different 
magnitude ranges, as listed in Table [T] 

At each epoch a new realization of the counts is made. Counts 
from stars, smooth backgrounds and instrumental backgrounds 
(thermal background and dark current) are Poisson realized, and 
fluctuations from readout noise are Gaussian realized. Photometry 
is performed on both the realization and a "true" image in a small, 
square, 3x3 pixel apertur^ around the microlensing source with 
the true (input) value of the smooth background subtracted, i.e., the 
number of counts from stellar sources in the aperture is measured 
to be 

JVpix 

iV= ^ (iVtot,.-{iVbg)), (6) 

i 

where Ntot,i is the total number of counts in pixel i and (A^bg) is 
the expectation of the counts in each pixel due to all the sources of 
smooth background, astrophysical and instrumental. An additional 
gaussian fluctuation of variance {asysN)^ is added to the final re- 
alization of the photometric measurement to simulate the effect of 
a systematic error floor. The photometric error is calculated as 

iVpix 

^N^^ {Ntot,i + O-rcad) + A^^f^sys, (7) 

i 

where A'^pix = 9 is the number of pixels in the aperture. The 
of the realized photometry relative to the "true" photometry is the 
of the true model which is used to calculate the Ax^ detec- 
tion statistic (see the next section). Should the number of counts 
in a pixel (including an assumed bias level) exceed the full well 
depth of the detector, then the pixel saturates. If that pixel lies in 
the photometry aperture the data point is removed from further cal- 
culations. 

It can be argued that the aperture photometry we simulate here 
is not appropriate for crowded fields, and that some form of PSF fit- 
ting photometry would be more realistic. While it is the case that 
the photometry that is performed on Euclid data will utilize the 
well known properties of the PSF to increase the photometric accu- 
racy, it should be noted that over the small number of pixels where 
we perform photometry, a boxcar is a reasonable approximation 
of the undersampled PSF. To check that the photometric method 
we use does not significantly impact the number of planet detec- 
tions we ran a test simulation performing photometry over a larger 
aperturelfl weighting the number of counts in each pixel by the in- 
tensity of the PSF in that pixel - this weighting approximates the 
performance of PSF fitting photometry. For a _ff-band survey and 
Earth-mass planets we find that weighted photometry increases the 
number of planet detections by 6 ± 2 percent. The improvement 
will be larger for less massive planets and smaller for more mas- 
sive planets, but in all cases will be too small to significantly affect 
our results. The improvement will be smaller for all other bands, 
because the smaller PSF in each case, and the smaller pixels on the 
VIS instrument, reduces the effect of blending. Indeed, the small 
under prediction of planet yields will likely be compensated by ef- 
fects that we do not model in our simulations (e.g. cosmic rays 
or common problems that affect infrared arrays such as ghosting. 



Testing .showed that the 3x3 aperture was the optimum aperture size for 
simple aperture photometry in our crowded fields. 

^ The radius of the aperture was 0.92", covering 29 pixels. This was chosen 
by experimentation to optimize the photometry. 
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NISPH NI5PJ NISPY VIS RIZ 

Figure 4. Top left: Example of a simulated false-colour composite image of a typical stai-field from the Euclid MAB/iLS simulation, with colours assigned 
as red-NISP H, green-NISP J and blue-VIS RIZ, each with a logarithmic stretch. The light green box surrounds the region that is shown zoomed-in in 
lower panels. The image covers 77 X 77 arcsec, equivalent to 1/64 of a single NISP detector, of which there are 16. These are shown to the right. Top 
right: Approximate representation of the NISP instrument 'paw-print'. The white areas show active detector regions, while black areas show the gaps between 
detectors. In the comer of one of the detectors is shown the size of a simulated image relative to the detectors. Bottom panels: The bottom panels show a small 
image region sun'ounding a microlensing event (located at the center and marked by cross-hairs), the top row showing images at baseline and the bottom row 
showing images at peak magnification ^ = 28. Panels from left to right show NISP H, J, Y, and VIS RIZ images, respectively. The small red box and red 
circle show the size of the aperture that was used to compute photometry in the NISP and VIS images respectively. 



charge diffusion or non-uniform pixel response functions), which 
are far more likely to degrade photometry than improve it. 

The properties of the detector/filter combinations that we 
model are listed in Table |2] We note the following about the pa- 
rameters listed in the table: 



• The zero-point is the AB magnitude of a point source, which 
would cause one count s^^ in the detector, after all telescope and 
instrument inefficiencies have been accounted for. The Euclid zero- 
points assume end-of-life instrument performance (M. Cropper, G. 
Seidel, private communication). 
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Table 2. Parameters of the Euclid telescope and detectors. Unle s s foot noted, 
all parameter values have been drawn from the iLaureiis et al.l i201 ih . Val- 
ues in brackets are values adopted for a longer slew time of 285 sec rather 
than our baseline assumption of 85 sec. Where necessary parameters are 
explained further in the text. 

Telescope parameters 

Diameter (m) 1.2 
Central blockage (m) 0.4 
Slew + settle time (s) 85(285) 



Detector parameters 



Instrument 


VIS 




NISP 




Filter 


RIZ 


Y 


J 


H 


Size (pixels) 


24k X 24k 




8k X 8k 




Pixel scale (arcsec) 


0.1 




0.3 




PSF FWHM (arcsec) 


0.18 


0.3* 


0.36* 


0.45* 


Bias level (e~) 


380t 




380t 




Full well depth (e~) 


216 




216 




Zero-point (ABmag) 


25.58* 


24.25** 


24.29** 


24.92** 


Readout noise (e~) 


4.5 


7.5* 


7.5* 


9.1* 


Thermal background 





0.26 


0.02 


0.02 


(e- s-i) 










Dark current (e~ s~^) 


0.00056* 




0.1* 




Systematic error 


o.ooit 




o.ooit 




Diffuse background 


21.5* 


21.3* 


21.3* 


21.4* 


(ABmag arcsec"^) 










Exposure time (s) 


540(270) 


90 


90 


54 


Images per stack 


1 


3(1) 


3(1) 


5(2) 


Readout time (s) 


< 85 




5* 





Schweitz er et al ]il2010l). The readout noise depends on the number of 
non-destructive reads; see text for further details. 

* Assumed in this work. 

*M. Cropper, private communication. 
**G. Seidel, private communication. 

*CCD203-82 data sheet, issue 2, 2007. e2v technologies, Elmsford, NY, 
USA. 

* Calculated base d on field locations, taking values for the zodiacal 
background from lLeinert et al 1l ll998h . and assuming an extt'a 0.2 
magnitudes from other sources such as scattered light. 

• We distinguish between dark current and thermal background. 
The dark current is the rate of counts induced by thermal sources 
within the detector pixels, and is independent of the observing 
band. The thermal background is the count rate due to thermal pho- 
tons emitted by all components of the spacecraft that hit the detec- 
tor, and is therefore affected by the choice of filter 

• For the Euclid simulations, we assume that the smooth back- 
ground is due primarily to zodiacal light. To account for any addi- 
tional smooth backgrounds we add an additional component with 
20 percent of the intensity of the zodiacal light. The zodiacal light 
background is calculated for e ach band at an elong ation of 90 ° in 
the ecliptic using data given bv lLeinert et al.l ( 119981) . 

• The VIS RIZ- and NISP F-bands are not included in the 
Besan§on model, so we assume that the AB magnitude of a star 
in the RIZ-hwiA is the average of its R and / AB magnitudes, and 
similarly we assume that the F-band magnitude is the average of / 
and J. 

Should a pixel within the photometry aperture saturate, the 
data point is flagged and is not included in the subsequent analysis. 
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1280 1285 1290 1295 1300 

Time/days 

Figure 5. Lightcurve of the simulated event shown in Figure |4] Fluxes are 
plotted normalized to the baseUne and blending in the ff-band. Grey, red, 
green and blue show data from NISP H, J, Y and VIS RIZ, respectively. 
The event reaches a peak magnification of ~28, but the normahzed flux 
only increases by a factor of ~3.3 because the source (-ffAB = 20.9) is 
blended with the diffraction spike of a much brighter star about 10 NISP- 
pixels away, and several other fainter stars, including the lens (Hab = 
32.0). At baseline the source contributes just 8 percent of the total flux in 
the H-band photometry aperture, though in the i?/Z-band aperture it con- 
tributes about 18 percent. Some of the event parameters are shown above 
the figure: Mi is the host-star mass; Ax^ is introduced in the next sec- 
tion. The inset shows the peak of the event, where a planetary signature 
is clearly detected, relative to the single-lens lightcurve (the grey line) that 
would have been observed were the planet not present. Data points are not 
scattered for claiity. 



We do not include the effects of cosmic rays in the images, except 
implicitly through the use of end-of-life instrument sensitivity val- 
ues. For the Euclid simulations, cosmic rays will only significantly 
affect observations with the VIS instrument, because the NISP in- 
strument, made up of infrar ed arrays, will use up-the-ramp fitting 
with non-destructive reads (lFixsenetal.ll2000l) to reduce readout 
noise and correct detector nonlinearities dSch weitzer et alj|2010l ; 
iBeletic et aLll2008l) . As a consequence of the multiple reads, up- 
the-ramp fitting mitigates against data loss due to cosmic rays and 
saturation. In order to ensure conservatism, we assume data with 
saturated pixels is lost completely. Currently we simulate the NISP 
instrument as a conventional CCD, but with variable read-noise de- 
termined by a fundamental read-noise (13 e~) and the number of 
non-destructive reads during an expos ure, which we a ssume occur 
at a constant rate of once every ~5 s dSchweitzer et al...2010) . We 
do not currently si mulate the more comp licated effects of charge 
smearing (see, e.g.. ICropper et al1l2010l) and ghosts from bright 
stars. 

For the Euclid simulations we use numerical PSFs computed 
for each instrument and each band. The NISP PSFs are computed 
near the edge of the detector field of view and include the effect of 
jitter and instrument optics in the worst case scenario (G. Seidel, 
private communication). The VIS PSF is similarly computed (M. 
Cropper, private communication). Figure |4]shows an example of a 
simulated, colour-composite image of a field with a microlensing 
event at its centre. The very brightest, reddest stars in the image are 
bright bulge giants of ~1 solar mass and ~80-120 solar radii. The 
much more numerous, but still bright and red, stars are red-clump 
giants in the bulge; bluer stars of a similar brightness are main- 
sequence F-stars in the disc. The fainter, resolved stars are turn-off 
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and upper-main-sequence stars in tiie bulge. Tiie figure also shows 
an approximate representation of the scale of the NISP instrument, 
which is constructed from 4x4 HgCdTe infrared arrays, each of 
2048 X 2048 pixels covering 10 x 10 arcmin, for a total detector 
area of 0.47 deg^ the gaps between detectors are approximately to 
scale. We do not include these gaps in the simulation and assume 
the instrument is a single 8192 x 8192-pixel detector. The lower 
section of Figure |4] shows a set of zoomed-in image sections, cen- 
tred on the microlensing event at peak and at baseline, in each of 
the NISP and VIS bands. Note the diffraction spikes and Airy rings 
in the VIS images; such spikes and rings can significantly affect 
photometry of faint sources. Figure |5] shows the lightcurve of the 
simulated event with peak magnification fi = 28 that occurs in the 
example image, including the points that are lost to saturation. For 
the sake of computational efficiency only a small image segment, 
just bigger than the largest aperture, is simulated in the standard 
operation of MAB/iLS. 

4.6 Planet detections 

To determine whether a bound planet is detected in a microlensing 
event we use three criteria, which will be further explained below: 

(i) the Ax^ between the best-fitting single-lens model and the 
best-fitting planetary model must be greater than 160, 

(ii) the Ax^ contribution from the primary observing band must 
be at least half of the total Ax^, 

(iii) the time of closest approach between the lens star and the 
source (to) must be within one of the 30-day observing seasons. 

For the first criterion, we assume the best-fitting planetary 
model to be the true underlying model that was used to simu- 
late the event. We choose Ax^ > 160, which corresponds to a 
a > 12.6 detection of the planet, because we find that the signa- 
tures of low-mass planets at this level of significance can usually 
be seen as systematic deviations from a single-lens lightcurve by 
eye (see e.g . even t (c) in Figure |7] below). This is in contrast to 
lGouldetalJ ( l2010h . who choose a higher threshold Ay^ > 500 
for pla nets in high-magnification microlensing events. iGould et"an 
were analyzing data from multiple, small ground based ob- 
servatories, which can suffer from various systematic effects (e.g. 
due to weather, differences in instrumentation, atmospheric effects 
in unfiltered data, etc.) that make the accurate estimation of photo- 
metric uncertai nties, and hence also extremely difficult. More 
recent work bv lYee et al.l ( l2012l) suggests that while a threshold of 
Ax^ > 500 may be appropriate for planets in high-magnification 
events, a lower threshold of Ax^ ^ 200 is likely to be more appro- 
priate for ground-based detection of planets in standard microlens- 
ing events, where the planetary signal is less ambiguous than in 
high-magnification events. A space-based microlensing data set 
will be much more uniform than ground-based data and will have 
much better characterized systematic effects, especially in the case 
of Euclid, whose design is driven by difficult, systematics-limited, 
weak lensing galaxy shape measurements. In order to see if plane- 
tary parameters could be measured from Ax^ ~ 160 lightcurves, 
we fitted a few simulated lightcurves using a Markov Chain Monte 
Carlo minimizer and found that even with Ax^ ~ 100 it was still 
possible to robustly measure the basic microlensing event param- 
eters, including the mass ratio, separation and in events where it 
was importa nt, the source radius crossing time (see Appendix A of 
|Pennvll201lh . 

Our choice of Ax^ > 160 also aids compa rison with other 
simulations which have chosen the same threshold jSennett & Rhid 




Figure 6. The effect of changing the Ax^ threshold on the number of 
planet detections. The number of planet detections with a Ax^ threshold 
^'^thrcsh' I'sl^tivs to 'he number of detections with Ax^ > 160, is plot- 
ted against Ax'^^^^^^^^■ Red, green, blue, magenta and cyan lines show the 
number of detections for 0.1-, 1-, 10-, 100- and 1000-A/g^ planets, respec- 
tively. 



I2OO2I Gaudi et al., unpublished), and is also the value adopted by 
the WFIRST science de finition team for th eir calculations of the ex- 
oplanet figure of merit jGreen et al.ll201 lb . Despite the widespread 
adoption of Ax'^ > 160 as a threshold for planet detections in 
space-based microlensing surveys, it is worth considering the ef- 
fect of changing the threshold. Figure |6] plots the number of de- 
tections with Ax^ greater than a variable threshold Axthrcsh> rel- 
ative to our chosen threshold of 160. A higher Ax^ threshold of 
Axthrosh — 200 would reduce yields by only ~25 percent for 
Mars-mass planets and less than ~10 percent for higher mass plan- 
ets. Even an e xtremely conservati ve threshold Axthrcsh = 500, 
such as used bv lGouldetaTI ( l2010l) for ground-based observations, 
reduces detections by 40-20 percent, depending on planet mass, 
above IA/q. Such a reduction in yield would not prevent Euclid 
from probing the abundance of Earth-mass planets, but may signif- 
icantly affect measurements for Mars-mass planets. However, such 
an extremely conservative cut will almost certainly not be neces- 
sary. 

Returning to the definition of selection criteria, the second cri- 
terion is chosen in order to allow fair comparisons between the dif- 
ferent bands that Euclid can observe in. By requiring that the con- 
tribution to Ax^ from the primary observing band is at least half 
of the total Ax^, we ensure that the primary band provides most 
of the information about the planet, and do not count as detections 
events where a planet is detected but most of the data are lost (e.g., 
due to saturation) or provides little information. 

The final criterion is chosen to increase the chance that the mi- 
crolensing event timescale is well constrained. The season length 
for microlensing observations on Euclid will be short, ~30 days, 
due to the restrictions of the spacecraft's sunshield. This can result 
in only a small portion of longer timescale events being monitored, 
and may also mean that the event timescale can not be constrained. 
To first order, it is the ratio of the timescale of the planetary per- 
turbation to the timescale of the main microlensing event which is 
used to measure the planetary mass ratio. Without the denominator 
of the ratio, the planetary mass ratio, and hence planetary mass can- 
not be estimated. Note however that it may be possible to constrain 
the event timescale from the curvature of the lightcurve without the 
peak, as in event (d) shown in Figure |7]below. 
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Figure 7. Example lightcurves from the MAB/iLS simulation of ExELS. The left column shows Earth-mass planet detections, with (a) showing a strong 
detection, (c) showing a detection close to the Ax^ threshold, and (e) showing an Earth-mass free-floating planet detection, (f) shows the full lightcurve of 
the free-floating planet detection in (e). The lightcurve in (b) shows a Mars-mass planet detection, but with the data points not scattered about the planetary 
lightcurve in order to emphasize the relative sizes of the photometric error bars. The lightcurve in (d) shows a O.OS-A/q planet that causes a signal well above 
our Ax^ threshold, but which is not counted as a detection because we require that the time of lens-source closest approach (the peak of the primary lensing 
event) be within an observing season. The inset figures, where included, either zoom in on planetary features or zoom out to show a larger section of the 
lightcurve. Grey, red, green and blue points with error bars show the simulated photometric data, while the black line shows the true lightcurve and the grey 
line shows the point-source single-lens lightcurve that would be seen if the planet were not present. In (e) the grey curve shows the lightcurve that would be 
seen if the source were a point, not a finite disc as is actually the case. In each hghtcurve the flux has been normaUzed to the i?-band flux, taking into account 
blending. The host mass, planet mass and semimajor axis, and Ax^ are shown above each lightcurve. 
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Figure |7] shows some example lightcurves from the simula- 
tion. The lightcurves show planet detections with varying degrees 
of significance, ranging from a detection that narrowly passed the 
Ax^ cut (lightcurve (c), Ax^ = 177) to a very significant detec- 
tion (lightcurve (a), Ax^ = 1527). Note however, that many events 
will have much higher Ax^ than this, up to Ax^ ~ 10®"^. The ex- 
ample lightcurves also cover a range of host and planet masses; the 
event with the lowest-mass planet is event (d), which has a planet 
mass Mp = O.O3i\f0 and detected with Ax^ = 384; however, 
due to our second criterion that to must lie in an observing season, 
this event is not counted as a detection. 



4.6.1 Free-floating planets 

To determine the expected number of free-floati ng planet de t ection s 
we adopt similar detection criteria to those of ISumi et al. I (l201lh . 
We require that in order to be classed as a detection, a free-floating 
planet lightcurve must have: 

(i) at least 6 consecutive data points (in any band) detected at 
greater than 3a" above baseline; and 

(ii) Ax'^ > 500 relative to a constant baseline model, using all 
the data points in the primary observing band that satisfy the first 
criteria. 

These criteria are in fa ct far more s tringe nt than the correspond- 
ing criteria imposed by ISumi et al] (l201lh . but we chose them to 
remain conservative, as we do not impose other criteria relati ng to 
the qu ality of microlensing model fits and images that iSumi et alj 
( I2OI II) use. 



5 EXPECTED YIELDS 

In this section we discuss the results of the MAB/^LS simulations of 
ExELS. Unless otherwise noted, we present the results assuming 
that each lens star in the simulation is orbited by a single planet of 
mass Mp with semimajor axis in the range 0.03 < a < 30 AU. 

Figure[8]shows the expected number of planet detections plot- 
ted against planet mass, using a naive assumption that there is one 
planet of mass Mp and semimajor axis 0.03 < a < 30 AU per star. 
The error bars on this plot, and all subsequent plots of the yield, 
show the uncertainty due to the finite number of events that we 
simulate. Error bars are not shown for the free-floating planet sim- 
ulation as they are similar to or smaller than the line thickness. For 
this naive assumption we expect a Euclid planetary microlensing 
survey would detect ~8, 38 and 147 bound Earth-, Neptune- and 
Saturn-mass planets (within 1-decade wide mass bins), and roughly 
half as many free-floating planets of the same masses. Euclid is sen- 
sitive to planets with masses as low as O.O3M0, but the detection 
rate for such low-mass planets is likely to be small unless the exo- 
planet mass function rises steeply in this mass regime. 

Recent measurements of planet abundances using several 
techniques have shown that the often used logarithmic planet 
mass function prior is quite unrealistic. Multiple studies have 
suggested t hat the number of planets increases with decreasing 
plane t mass jCurmning et alj2008l; Johnson et al .l20ld;l ISumi et aU 
l20ld ; [Howard et alj|2012l ; iMavor et al.ll201ll; ICassan et alj|2012l) 
and that planets are not distributed logarithmically in semimajor 
axis (Gumming et al. 2008). T his picture is also suppor te d by planet 
popul ation synthesis models l lMordasini et alj|2009al lbl Tlda & LinI 
I2OO8I) . In Figure |9] we consider a more realistic two-parameter 
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Figure 8. Number of planets detected in a 300-day sureey by Euclid, plotted 
against planet mass Mp. The survey is primarily conducted in the NISP H- 
band. The solid black line shows the expected bound planet yield, assuming 
one planet of mass Mp per star with semimajor axis 0.03 < a < 30 AU; 
error bars show our estimated statistical errors from simulations of a finite 
number of Ughtcurves. The solid grey line shows the yield if the third cut on 
the time of the event peak is not appHed. The dashed line shows the expected 
yield of free-floating planets, assuming there is 1 free-floating planet per 
Galactic star. The masses of Solar System planets are indicated by letters, 
and the numbers above/below the lines show the yields when applying the 
full sets detection criteria. 



power-law planetary mass function: 

f(M ) = '^'^ - f (Mz 
^^''^P^ - dlogMpdloga '^'{m. 



(8) 



where f{Mp) is now the number of planets of mass Mp per decade 
of planet mass per decade of semimajor axis per star and where /. 
is the planet abundance (in dex~^ star~^) at some mass M, about 
which the mass function pivots. Here, a is the slope of the mass 
function, with negative values implying increasing planetary abun- 
dance with decreasing planetary mass. For simplicity, and because 
there are no measurements of the slope of the planetary semimajor 
axis distributions in the regime probed by microlensing, we assume 
that dN/d log a is constant. 

We use two estimates of the mass-function parameters based 
on measurements made using both RV and microlensing data sets. 
The first, more conservative mass function (in terms of the yield of 
low-mass pla nets) uses the ma ss-function slope a = —0.31 ± 0.20 
measured bv ICumming et al" I Jioos) from planets with periods in 
the range T = 2-2000 d, detected via radial velocities. For the 
normaliz ation we use /. = 0.36 ± 0.15 at M, ~ 8OM0, mea- 
sured by iGould et al.l ( 120 lOl) from hifih-magnifj cation microlens- 
ing events observed by MicroFUN. lOould et"al] ( I2OIO ') argue that 
this value is consistent with the abundance and semimajor axis dis- 
tribution measured by ICumming et al. I ( l2008h . extrapolated to or- 
bits with a ~ 2.5 A U. We note that the host stars studied by 
ICumming et"ai] ( |2008|) typically have higher masses than those 
that are probed by mi crolensing. We call the combin ation of the 
ICumming elal] ( 120081) slope and iGould et al] ilOld) normaliza- 
tion, the RV mass function. The second mass function we consider 
has a mass function slope a = —0.73 ± 0.17 and norma lization 
/. = 0.24l°;ig at M, = 95M0, as measured by Cas san et alj 
( I2OI2I) from microlensing detections. We call this the microlensing 
(/iL) mass function. We note that at low masses, the extrapolation 
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Table 3. Expected total number of planet detections by a 300-day Euclid 
microlensing survey for different mass functions (with planet masses in the 
range 0.03 < Mp/Mgj < 3000 (roughly 0.6 Mercury-mass to 10 Jupiter- 
mass). 



Mass function 


Number of detections 


log-log 


718 


RV 


502 




541 


/iL saturated 


356 



of the microlensing mass function implies close packing of plane- 
tary systems. We also plot the microlensing mass function assum- 
ing that it saturates at a planet abundance of 2 dex~^ star~^. How- 
ever, we note that the Keple r 20 planetary syste m comprises five 
exoplanet candidates so far ( lGautieretal1l2012h . all within about 
1 dex in both mass and separation. Our saturation limit is therefore 
likely to be conservative. 

Figure |9] plots the yields that would be expected for the dif- 
ferent mass functions. Perhaps the most important thing that the 
top panel of Figure|9]highlights is that, despite the degree of uncer- 
tainty in the extrapolation to low planet masses provided by empiri- 
cal estimates of the mass functions, we can expect a 300-day Euclid 
survey to detect significant numbers of planets of Mars-mass and 
above. Table [3] shows the total number of detections expected for 
each mass function. The number of expected detections imply that 
Euclid data would allow the different model mass functions to be 
discriminated between. In fact, we can look at the power of Eu- 
clid to measure the mass function more easily in the lower panel of 
Figure |9] 

The lower panel of Figure|9]shows the expected uncertainty on 
the planet abundance in one-decade mass bins, assuming the satu- 
rated microlensing mass function and that half the Euclid planet 
detections have mass measurements. Such mass measurements can 
be made by estimating the mass of the host from photometry of 
the host star. Such an estimate should be possible for many of the 
hosts using ExELS survey data alone, thanks to the extremely deep, 
high-resolution images that can be built by combining the randomly 
dithered survey images. A stack of such images, one for each sea- 
son, will allow the light of the source, host and any unrelated stars 
to be disentangled as they separate due to their mutual proper mo- 
tions. [Bennetretal] i2007h give a detailed discussion of how these 
mass measurements are made. ISennett et al.l ( l2007h estimate that 
such mass measurements should be possible in most space-based 
planetary microlensing detections. However, it may be the case 
that the larger pixels of the NISP instrument preclude full photo- 
metric host-mass measurements, but even if this is so, the deep, 
high-resolution images from the VIS channel should provide con- 
straints. Even without mass measurements. Figure |9] indicates the 
precision of measurements of the mass-rafto function, which would 
encode much of the same information. The uncertainties shown by 
the coloured bands in the figure are the uncertainty on the absolute 
abundance of planets in each bin. This i s in contrast to measure- 
ments such as those of Gum ming et alj (2008) and Cassan et al. 
( |2012|) . which are the uncertainties on a small number of power- 
law model parameters assuming the models are correct. If we here 
assume that the saturated microlensing mass function is correct, 
then we can see that a 300-day Euclid microlensing survey would 
measure the abundance of Earth mass planets to be 2 per star with 
a significance of 4.7-(7, and similarly the abundance of Mars mass 
planets to 1.7 -a. However, if the microlensing program were only 
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Figure 9. The upper panel shows predictions of the planet yield based on 
recent estimates of the planet abundance and planet-mass distribution. The 
solid line shows a naive logarithmic prior of one planet per decade of mass 
and semimajor axis per star. The dashed line (labelled RV) shows the ex- 
pected yield using an e xtrapolation of the mass-function slope measured by 
ICumming etaP i200^ using RV data combined with a normaUzation mea- 



from microlensing data, which iGould et alj 



rig 

sured by IGould et al. _ 

iioiol) argue i s compatible wi t h the slope of the semimajor axis distribu- 
tion found by ICumming et al.l )2008l) . The dot-dashed line (labelled /xL) 
shows the expected yield using the mass function slope and normalization 
measured from microlensing data bv lCassan et al ]l|2012). A branc hing dot- 
dashed line, and the numbers above it, show the yield if the Cass an et alj 
(2012) microlensing mass function saturates at 2 planets per dex'^ per star. 
The lower panel shows the actual forni of each of the mass functions shown 
in the top panel. The filled, coloured regions show the size of model- 
independent 1-(T statistical (square root A'^) errors on measurements of the 
planet abundance in 1-decade mass bins centred at Mp, assuming the sat- 
urated microlensing mass function and also assuming that only half of the 
planet detections have host mass measurements. The red, green, magenta, 
cyan and grey regions show the error bars for Euclid microlensing surveys 
lasting 60, 120, 180, 240 and 300 days respectively. This implies a 300- 
day Euclid microlensing survey would measure the abundance of Earth- 
mass and Mars-mass planets to approximately 4.7 and 1.7-(t respectively, 
whereas a 120-day Euclid survey would reach just 3.0 and 1.2-ct signifi- 
cance, both assuming the saturated mass function. 
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Figure 10. The sensitivity of Euclid in tlie A/p-a plane. Red lines show the 
expected yield of a 300-day Euclid sui'vey with 60 days of observations per 
year, plotted against planet mass and semimajor axis, assuming one planet 
per star at each point in the planet mass-semimajor axis plane. Horizontal 
arrows are plotted when the expected yield of free floating planets of that 
mass exceeds the yield of bound planets (assuming one free floating planet 
per star). The grey points show planets listed by the Exoplanets Orbits 
Database as of 17th March 2012 ( Wright et al.l l201 iV and Ught blue points 
show candidate planets from the Kepler mission iBatalha et al. 2013), with 
masses calculated using the mass-radius relation of Lissauer et al.l l l201 II) . 
The red points show planets detected via microlensing to date. 

120 days, the significance of the abundance measurements would 
reduce to 3.0- and 1.2-cr for Earth- and Mars-mass planets, respec- 
tively. 

5.1 The Mp-a diagram 

We have discussed the ability of our simulated survey to probe 
the planetary mass function, but a perhaps more important goal of 
such a survey is to explore the planet mass-semimajor axis (Mp-a) 
plane whe re planet f ormation models pre dict a lot of structure (e.g., 
Ilda & Lin 2004; Mo rdasini et al.ll2009al) . Figure [To] plots contours 
of planet detection yields for the simulated survey in the Mp-a 
plane, assuming there is one planet per host at a given point in 
the plane. The positions of planet detections to date, by all detec- 
tion methods (RV, transits, direct detection, timing and microlens- 
ing) are also shown, as well as candidate planets detected by Ke- 
pler teatalha et alj2013h , which have been plotted by assuming the 
planetar y mass-radius relation , Mp = (i?p/i?0)^'''^M0, which is 
used bv lLissauer et al ll). It is clear that microlensing surveys 
probe a different region of the Mp-a plane to all other detection 
methods, covering planets in orbits ~0.2~20 AU, as well as free 
floating planets. Note that microlensing can be used to detect plan- 
ets with any semimajor axis larger than ~20 AU, but there is a sig- 
nificant chance that the microlensing event of the host will not be 
detectable. These cases will be classified as free-floating planet de- 
tections (see e.g. jSumi et al. 2011; Bennett et aL 2012). The peak 
sensitivity of the simulated Euclid survey is at a semimajor axis 
a »i 1.5-5 AU, in good agreeme nt with previous simu lations of 
space-based microlensing surveys jBennett & Rhiell2002l Gaudi et 
al., unpublished). The planets to which Euclid is sensitive lie in 
wider orbits than those detectable by Kepler, and stretch to much 
lower masses than can be detected by RV in this semimajor axis 
range, reaching down to Mars mass. The range of semimajor axis 
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Figure 11. Predictions of the planet yield as a function of semi-major axis a. 
The red, green, blue, magenta and cyan lines denote yields for 0.1,1,10,100 
and 1000 Mgj, respectively. 
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Figure 12. The nuinber of planet detections plotted against the planetary 
effective temperature, assuming an albedo of 0.3. Lines are as for Figure lTTI 

probed by Euclid decreases with decreasing mass, from ~0.2 to 
more than 20 AU for Jupiter-mass planets, down to ~1-14 AU for 
Earth-mass planets and ~1.5-5 AU for Mars-mass planets. There 
will be a significant degree of overlap between Euclid and full- 
mission Kepler detections at separations 0.3 a < 1 AU. Sim- 
ilarly, at masses larger than Mp > 5OA/0, there will be overlap 
with RV surveys over a wide range of semimajor axes. Both over- 
laps will facilitate comparisons between the data sets of each tech- 
nique. It should be noted however, that the host populations probed 
by each technique are different, as we will see in the next section. 

Figure [TT] plots the expected yield for various planet masses 
as a function of semimajor axis a, using a simplistic assumption 
of one planet per host at the given mass and separation. The peak 
sensitivity of Euclid is to planets with semimajor axis a « 1.5- 
5 AU. The sensitivity is ~10 percent of the peak sensitivity in the 
range 0.5 < a < 20 AU. 

Figure [T2]plots the distribution of planet detections as a func- 
tion of the effective temperature of the planet, calculated as 

Tcff,p = yj(i-A)^/*reff,i, (9) 

where Ri is the radius of the host star, A is the planet's albedo, as- 
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Figure 13. The distribution of microlensing timescales. Red curves sliow 
the distribution of event timescales for stellar microlensing events. The 
solid curve shows all events with impact parameter uq < 1, regardless 
of whether they are detected, the dashed curve shows events which are de- 
tected above baseline with Ax^ > 500, and the dot-dashed curve shows 
those detected events which peak during an observing season. The solid 
grey hues show the theoretically expect ed asymptotic slope of the distri- 
bution, with power law slopes of ±3 ( Mao & Paczvnsld[l99^ . The cyan 
data points show the timescale distribution observed by the MOA sur- 
vey jSumi et al]|201lh , which is uncorrected for detection efficiency and 
scaled arbitrarily - this is most closely comparable to the dashed line show- 
ing all events detected by Euclid. The black lines and data points show the 
timescale distribution for events with detected planets. The solid line shows 
the timescale distribution of the host star microlensing event for 100- A/q 
planet detections with no restriction on to, while the data points show the 
same, but only for events where to lies in an observing season. The dashed 
and dot-dashed lines show the timescale distribution of detected IOO-Mq 
and free-floating planet detections, respectively. 



sumed to be = 0.3 and Tcs,i is the effective temperature of the 
star. Both Ri and Teff,i are provided as outputs of the Besanjon. 
The distribution of detected planet temperatures peaks at ~ SO- 
SO K, with a long tail towards lower temperatures and a rapid de- 
chne towards higher temperatures. However, there should still be a 
small number of detections of planets with effective temperatures 
> 200 K. 



5.2 Host star properties 

The primary observable of the microlensing lightcurve that is re- 
lated to the host star's mass is the event timescale. The timescale is 
a degenerate combination of the total lens mass, the relative lens- 
source proper motion and the distances to the source and lens. 
Figure [13] plots the timescale distributions of all the microlens- 
ing events that occur within the observed fields and also the dis- 
tributions for several cases of planet detections. The timescale dis- 
tribution for bound planet detections is similar to the underlying 
timescale distribution, but is affected by the choice of detection cri- 
teria. Our third criterion, designed to select only events with well 
characterized timescales, cuts out potential detections in some long 
timescale events. Some of these events are detections of planets 
with large orbits, where the planetary lensing event is seen but the 
stellar host microlensing event is only partially covered (in which 
case the planet parameters may be poorly constrained) or may be 
missed completely (in which case the planet event would enter the 
free-floating planet sample). However, in other cases the cut on to 
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Figure 14. Predictions of the 100-A/gj planet yield as a function of lens 
(solid lines) and source (dashed line) distances, Di and Db, respectively. 
The red and green lines show the contributions due to bulge and thin disc 
lenses, respectively; thick disc and halo lenses contribute the remainder, 
which is small. 



is too zealous, and long timescale events with to outside the observ- 
ing window, but with significant magnification in several seasons, 
are cut from the sample. The timescale of these events, and hence 
also the planetary parameters, are likely to be well constrained. 

Figure[T3]also plots the free-floating planet timescale distribu- 
tions for planets of 1 and 100 Earth masses. Free-floating planets 
will dominate the timescale distribution at timescales less than a 
few days, if they exist in numbers similar to those suggested by 
ISumi et al.l ( l201ll) , which is twice the abundance that we have as- 
sumed. 

Figure [T4]plots the distribution of IOO-Mq planet detections 
as a function of lens and source distances, Di and Dg, respectively. 
The contribution of thin-disc and bulge populations to the yields is 
also plotted. Thick disc and stellar halo lens yields have not been 
plotted as at no point are they dominant. However, near the Galactic 
centre it should be noted that stellar halo lenses have a higher yield 
than the thin disc due to the disc hole (see Section l4Tt . Most of 
the host stars are near-side bulge stars between 5.5 < Di < 8 kpc. 
Beyond this, the number of lenses with detected planets drops-off 
exponentially with increasing distance, dropping by four orders of 
magnitude from 73i ~ 9 to 15 kpc. The steepness of this fall is 
partly caused by the truncation of the source distribution at 15 kpc. 
Though the majority of lenses are in the bulge, a substantial number 
reside in the near disc. The contribution of planet detections by each 
component is 68, 27, 1.2 and 3.5 percent for the bulge, thin disc, 
thick disc and stellar halo populations, respectively. Unlike the lens 
stars, the majority of source stars reside in the far bulge, with a 
small fraction in the far disc. Very few near disc stars act as sources 
due to the low optical depth to sources on the near side of the bulge. 

Figure [T5]plots the distribution of IOO-Mq planet detections 
as a function of the host star spectral type. The majority of hosts are 
M dwarfs, but there are a significant number of detections around 
G and K dwarfs and also white dwarfs. There will be a negligi- 
ble number of detections around F and earlier-type stars due to 
their low number density. The distribution of planetary host stars 
probed by Euclid is very different to that probed by any other tech- 
nique. For example, most of Euclid's host stars are M dwarfs in the 
bulge , whereas most of K epler's host stars are FGK dwarfs in the 
disc ( iHoward et alj201^ . 
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Figure 15. Histogram of the number of lOO-Afgj planet detections plotted 
against tlie effective temperature of the host star, binned according to the 
spectral type designations in the Besanfon model. The shaded region shows 
the contribution due to main-sequence host stars, while white regions show 
the contiibution of evolved host stars. The dotted line shows t he distribution 
of eff ective temperatures of high-priority Kepler target stars teatalha et alj 

l2oToh . 



6 VARIATIONS ON THE FIDUCIAL SIMULATIONS 

In the previous section we have investigated the potential planet 
yield of a Euclid microlensing survey and the properties of de- 
tectable planets and their hosts. In this section we investigate how 
the planet yield is affected by our choice of primary observing 
band, the level of systematic photometry errors and the choice of 
spacecraft design. 



6.1 Primary observing band 

We begin by examining the choice of primary observing band. The 
survey strategy we have simulated involves the majority of observa- 
tions being taken in a primary band with a cadence of ~18 minutes 
while auxiliary observations to gain colour information are taken 
every ~12 hours. We consider the use of each band available to 
Euclid, Y, J and H in the near infrared using NISP and RIZ using 
VIS. As NISP and VIS can image the same field concurrently we 
also consider simultaneous observations in RIZ and H. To main- 
tain a comparable cadence, when RIZ is the primary band (or VIS 
is operating simultaneously with NISP), the VIS exposure times are 
270 s, as opposed to 540 s when RIZ is used as an auxiliary band. 
In each scenario the total exposure time is identical, but the actual 
cadence is slightly different due to differences in the number of 
stacked images (we assume a 5 s overhead between the images in 
the NISP stacks, and the shutter on VIS takes 10 s to open or close). 
As the sensitivities of the instruments in each band are slightly dif- 
ferent, the images have different depths. 

Figure [16] shows the expected planet yields as a function of 
the primary observing band. Focussing first on the scenarios with- 
out simultaneous imaging, it is clear that _ff-band offers the high- 
est planet yields compared to the other two infrared bands. This 
is partly due to the increased depth from a stack of 5 images for 
H as opposed to a stack of 3 images for J and Y (the individual 
exposur e times have been ch osen optimize Euclid's cosmological 
surveys; iLaureiis et al.ll20ll1) . However, it is also due to the lower 
extinction suffered in the _ff-band, and the correspondingly higher 
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Figure 16. Expected planet detections plotted against the different primary 
observing bands. Free-floating planet detections are not included. 



number density of sources with magnitudes lower than the source 
catalogue cut-off of JJvcga < 24. 

The survey imaging with both available instruments simulta- 
neously obviously performs better than when using each instrument 
on its own. The increase in yield is ~22±4 percent for both Saturn- 
mass and Earth-mass planets. As for the single primary instrument 
scenarios, we require that the Ax^ contribution of the primary 
bands (the sum of RIZ and H) to be greater than half the total 
Ax^. In reality, the expected yield of the simultaneous imaging 
scenario represents an upper limit, as there are a number of lim- 
itations that may preclude simultaneous imaging with VIS for all 
pointings. These include losses due to cosmic rays, wh ich will af- 
fect ~20 percent of VIS data points (Laureii s et al .l201l l). downlink 
bandwidth and power consumption limitations, which may only al- 
low a simultaneous VIS exposure every other pointing, say. The 
increase in yield may therefore be small. However, the real value 
of simultaneous VIS imaging will be the increased number of ex- 
posures it is possible to stack in order to detect the lens stars. This 
will greatly increase the depth of VIS images stacked over the entire 
season, which in turn will allow the direct detection of more lens 
stars, and hence an increase in the accuracy and number of mass 
measurements it is possible to make. We discuss this further in Sec- 
tion|7] Simultaneous VIS imaging will also allow source colours to 
be measured in many low-mass free-floating planet events, which 
will help to constrain their mass. It is clear therefore that as many 
simultaneous VIS exposures should be taken as possible. 



6.2 Systematic errors 

There are many possible sources of systematic error, which can in- 
clude image reduction, photometry, image persistence in the detec- 
tor, scattered light, temperature changes in the telescope and source 
and intrinsic variability in the source, lens or a blended star. The 
magnitude and behaviour of each systematic will also be differ- 
ent; for example, temperature changes will likely induce long-term 
trends in the photometry, while image persistence may introduce 
a small point-to-point scatter together with occasional, randomly- 
timed outliers. It is likely that the systematics that produce long 
term trends may be corrected for, to a large extent, either by using 
additional spacecraft telemetry or by detrending similar to that used 
in transiting exoplanet analyses (e.g., lHolman et alj201(]|) . Even for 
some systematics that behave more randomly, it may be possible to 
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Figure 17. Expected planet detections plotted against the size of the sys- 
tematic ertor component. 
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Figure 18. Expected planet detections plotted for 85-sec and 285-sec slew- 
ing times, which encompass the likely range anticipated by different designs 
for the Euclid manoeuvring system. 



account for and correct errors; for example, it may be possible to 
correct for image persistence errors to a certain degree by u sing pre- 
ceding images. It is therefore difficult to predict the magnitude and 
behaviour of systematics a priori. We therefore choose to model 
systematic eiTors by assuming them to be Gaussian, and add the 
systematic component in quadrature to the standard photometric 
error While likely a poor model for the actual systematics, it ef- 
fectively introduces a floor below which it is not possible improve 
photometry by collecting more photons. 

Figure[T2]plots the expected planet yield against differing val- 
ues of the systematic error component that we assume. In all other 
simulations we have used the fiducial value of the fractional sys- 
tematic error fjsys = 0.001. Reducing the systematic error further 
from this point does not provide a significant increase in yield, as 
for the most part, at this level of systematic, photometric accuracy is 
limited by photon noise. Increasing the systematic to ctsys = 0.003 
does cause a drop in yields, by ~8 percent for giant planets to 
~20 percent for low-mass planets, as the systematic component 
becomes comparable to the photon noise. The situation is worse 
still for (Tsys = 0.005, where the systematic component dominates. 
However, even with a systematic error component this large and 
the conservative log-log mass function, ~4-5 Earth-mass planet 
detections can be expected. 

It is not possible at this stage to estimate the magnitude of sys- 
tematic error that should be used in our simulations, but it should 
be noted that ground-based microlensing analyses often have sys- 
tematic errors of a similar magnitude to the values that we have 
simulated (D. Bennett, private communication). The tight control 
of systematics required by Euclid for galaxy-shape measurements 
should mean that Euclid-VIS will be one of the be st-characterized 
optical instruments ever built dLaureiis et al.ll20l"ll) : similarly, NISP 
will be optimized for performing accurate, photometry of faint 
galaxies. Furthermore, Clanton et al. (2012) recently showed that 
the HgCdTe detectors that will be used in NISP can perform sta- 
ble photometry to ~50 parts per million. How these considerations 
will relate to crowded-field photometry is not yet clear, but it is al- 
most certain that the systematics will be lower than those achieved 
from the ground, potentially by a large factor Our fiducial choice 
of a fractional systematic error 0.001 (1000 parts per million) is 
therefore almost certainly conservative. 



6.3 Slewing time 

Another uncertainty in the yields we predict results from uncertain- 
ties in the spacecraft design. 

Whilst the manufacturer and final design for the Euclid space- 
craft is yet to be decided it is possible to explore some factors which 
are likely to have an important bearing on its microlensing survey 
capabilities. One important factor is the choice of manoeuvring sys- 
tem used for slewing between fields. For fixed exposure and areal 
coverage the slew and settle time determines the cadence it is possi- 
ble to achieve on a particular field. Alternatively for larger slewing 
times one may shorten the exposure time to maintain cadence and 
areal coverage. Since the detection of low mass planets depends 
crucially on cadence this alternative approach is preferable when 
considering the impact of adjustments to the slewing time. 

The slew time will ultimately depend on the technology used, 
in particular whether gas thrusters or reaction wheels are used to 
perform field-to-field slews. A plausible range for the slew times 
based on initial design proposals is 85-285 sec. 

Figure [18] shows the expected yield at either end of this slew 
time range. Maintaining a constant cadence of around 20 mins be- 
tween repeat visits to a given field allows 270 s per pointing of 
stacked //-band exposure time for 85-sec slews and 108 s for 285- 
sec slews. The increased depth allowed by a shorter slewing time 
produces a yield that is higher by 50 ± 12 percent at Earth mass 
and 22 ± 5 percent at Saturn mass. 



7 SUMMARY DISCUSSION 

The Euclid dark energy survey, which has been selected by ESA 
to fly in 2019, is likely to undertake additional legacy science pro- 
grams. The design requirements of the Euclid weak lensing pro- 
gramme also make it very well suited to an exoplanet survey using 
microlensing and the Euclid Exoplanet Science Working Group has 
been set up to study this proposal. 

We have developed a baseline design for the Exoplanets 
Euclid Legacy Survey (ExELS) using a detailed simulation of 
microlensing. The simulator, dubbed MABuLS, is based on the 
Besanjon Galaxy model (Robin etal. 20()3[). It is the first mi- 
crolensing simulator to generate blending and event parameter dis- 
tributions in a self-consistent manner and it enables realistic com- 
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parisons of the performance of Euclid in different optical and in- 
frared pass-bands. We hiave used MABELS to study a design for 
ExELS with a total observing baseline of up to 300 days and a sur- 
vey area of 1.6 deg^. We show that of the band-passes available 
to Euclid a survey primarily conducted in H will yield the largest 
number of planet detections, with around 45 Earth-mass planets 
and even ~6 Mars-mass planets. These numbers are based on cur- 
rent extrapolations of the exoplanet abundance determined by mi- 
crolensing and radial velocity surveys. Such low-mass planets in 
the orbits probed by Euclid (all separations larger than ~ 1 AU) 
are presently inaccessible to any other planet detection technique, 
including microlensing surveys from the ground. 

While space-based microlensing offers significantly higher 
yields per unit time than do ground-based observations, this is 
not the only motivation for space-based observations. A standard 
planetary microlensing event does not automatically imply a mea- 
surement of planet mass or semimajor axis, only the planet-star 
mass ratio and the projected star-planet separation in units of the 
Einstein radius te. To measure the planet mass we must mea- 
sure the lens mass, either by detecting subtle, higher-order ef- 
fects in t he microlen sing lightcu rve, s uch as microlensing paral- 
lax (e.g. iGould [200 0: An et all l2002h. or directl v detecting the 
lens star ('Alcock et al. 2001l; iKozlowski et alj2007l) . Wit hout these 
the mass can only be deter mined probabilistically (e.g., iDominilJ 
12004 iBeaulieu et al.|[20o3) . The projected separation in physical 
units can be determined if the lens mass and distance are known 
(as well as the source distance, which it is possible to estimate 
from its colour and magnitude). Determinin g the semimajor axis 
will require the dete ction of orbital motion dSennett et al.ll2010hl : 
ISkowron et al. 2011), but this will only be possible in a subset of 
events (IPennv etalj201lh . For a survey by Euclid we expect paral- 
lax measurements to be rare. Parallax effects are strongest in long 
microlensing events lasting a substantial fr action of a year due to 
the acceleration of the Earth (lGouldlll992h . but Euclid's seasons 
will be too short to cons train or detect a parallax signal in most 
events temith et alfcoOSh . 

However, thanks to the high-resolution imaging capabili- 
ties of the VIS ins trument, lens detection should be common 
teermett et alj|200^ . In events where the light of the lens is de- 
tected, the lens mass and distance can be determined by combin- 
ing measurements of the angular Einstein radius 9^ (which gives a 
mass-distance relation) with a main-sequence mass-luminosity re- 
lation. Measurement of 9^ should be possible for a large share of 
detected events, either from finite-source effects in the lightcurve 
or by mea suring the relative le ns-source proper motion as the pair 
separates (^Bennett et alj|20()7h . It is also possible to estimate the 
lens m ass and distance from measurements of its colour and mag- 
nitude teennett et"^l2007h . From a single epoch of NISP and VIS 
images, this will likely not be possible. However, over each 30-day 
observing period around 2000 images will be taken in NISP H- 
band, with possibly a similar number with the VIS camera. These 
images will have random pixel dither offsets. The images can there- 
fore be stacked to form a much deeper, higher-resolution image in 
each band. From these images it should be possible to isolate the 
source (whose brightness is known from the lightcurve) from any 
blended light. After subtracting the source, if the remaining light 
is due to the lens, its mass can be estimated from its colour and 
magnitude. The planet mass can then be determined, as the planet- 
host mass ratio is known from the lightcurve. However, if either the 
source or lens has a lu minous companion, estimating the lens mass 
will be more difficult ( iBennett et al.ll2007l) . 

We have not attempted to estimate the number of planet detec- 



tions with mass measurements in the present work, but we aim to 
study this in a future work. These calculations will allow a full de- 
termination of planetaiy microlensing figures of merit, such as the 
one d efined by the WFIRST Science Definition Team jOreen et alj 
[20Tlh . 

Finally, it is worth stating that our simulation of ExELS has 
not been optimized. There are many factors that can be varied to in- 
crease planet yields, such as the choice of target fields, the number 
of target fields and the strategy with which they are observed. How- 
ever, planet yields are not the only measure of the scientific yield 
of the survey. For example, planetary-mass measurements without 
the need for additional follow-up observations would be an impor- 
tant goal of the Euclid microlensing survey, and so any assessment 
of the relative performance of different possible surveys must also 
evaluate performances in this respect. 

We have shown that ExELS will be unrivalled in terms of 
its sensitivity to the cold exoplanet regime. A survey of at least 
six months total duration should be able to measure the exoplanet 
distribution function down to Earth mass over all host separations 
above I AU. This will fill in a major incompleteness in the current 
exoplanet discovery space which is vital for informing planet for- 
mation theories. This together with ExELS's ability to detect hot 
exoplanets and sub-stellar objects (Paper II) make it a very attrac- 
tive addition to Euclid's science capability. 
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